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NOMENCLATURE 


[A], 

= 

defined  in  Eq.  (19) 

[B1J 

= 

defined  in  Eq.  (19) 

[Cl. 

= 

defined  in  Eq.  (19) 

<C>j 

= 

defined  by  Eq.  (15) 

IF] 

= 

excitation  force  levels 

V 

= 

jl  element  of  [F] 

[G]] 

= 

defined  in  Eq.  (19) 

<G>] 

3 

jth  column  of  [G] 

3 

defined  in  Eq.  (19) 

(i>3 

3 

defined  by  Eq.  (12) 

L 

* 

Lagrange  function 

[M] 

= 

mass  matrix 

C 

3 

critical  damping  ratio 

[W] 

3 

weighting  matrix 

a.  . 
ij 

3 

defined  by  Eq.  (4) 

[B] 

3 

T 

eigenvalues  of  [$]  [M] [$] 

[9] 

3 

T 

eigenvectors  of  [$]  [M] [$] 

[X] 

3 

matrix  of  Lagrange  multipliers 

3 

element  of  Hi  of  [X],  Lagrange  multiplier 

'*>] 

3 

jth  row  of  [X] 
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NOMENCLATURE  (continued) 


[*] 

= 

all  the  normal  modes  of  a  structure 

[<H 

= 

identified  mode  shapes 

[?] 

= 

analytically  orthogonal i zed  mode  shapes 

[<l>ra] 

= 

measured  modes 

(♦}j 

= 

jth  column  of  [$] 

♦jk 

= 

jk  element  of  [$3 

A1" 

♦jk 

= 

jk  element  of  [$m] 

W 

= 

defined  by  Eq.  (17) 

[Q] 

* 

admittance  matrix,  imaginary  components 

Qii 

= 

ij  element  of  [fl] 

A 

[fl,i 

= 

defined  in  Eq.  (9) 

(l). 

J 

a 

frequency  of  excitation  (rad/sec) 

CO  . 

m 

= 

ith  natural  frequency  (rad/sec) 

0 

s 

element-by-element  matrix  multiplication  operator 
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I.  INTRODUCTION 


Accurate  structural  dynamic  models  of  complex  spacecraft  are  a  require¬ 
ment.  Unfortunately,  analytical  models  agree  closely  with  properly  measured 
mode  data  only  in  the  first  few  modes.  To  minimize  the  effects  of  this  defi¬ 
ciency,  two  approaches  are  widely  used.  The  first  approach  is  to  adjust  the 
analytical  dynamic  model  to  improve  correlation  between  the  analytical  and 
empirical  modes.  ,  Any  remaining  difference  between  the  two  mode  sets  is  then 
usually  ignored,  and  the  adjusted  analytical  model  is  adopted  as  the  model  of 
the  actual  hardware. 

The  second  approach  also  involves  adjusting  the  analytical  model  to  improve 
correlation  with  the  measured  modes.  However,  unlike  the  first  approach,  the 
measured  modes  are  then  adopted  as  the  normal  coordinates  of  the  dynamic  model. 
This  requires  that  the  mode  survey  test  article  be  representative  of  flight 
hardware.  The  principal  advantage  of  this  approach  is  that  the  influence  of 
deficiencies  remaining  in  the  analytical  model,  after  all  adjustments  have 
been  made,  are  minimized. 

Structural  dynamic  models  of  complex  space  systems  are  typically  formulated 
by  component  mode  synthesis  coupling  of  substructure  modal  models.  Inherent  in 
the  coupling  procedures  is  the  assumption  that  the  modes  of  a  substructure  are 
orthogonal  with  respect  to  its  mass  matrix.  Therefore,  if  measured  modes  are 
used,  any  deviation  from  mass-weighted  orthogonality  must  be  corrected. 

Numerous  procedures  have  been  proposed  to  analytically  orthogonalize 
measured  modes.  Gravitz  (Ref.  1)  proposed  calculating  a  symmetric  influence 
coefficient  matrix  by  averaging  the  off-diagonal  terms  of  a  matrix  obtained 
from  the  measured  frequencies,  measured  modes,  and  the  generalized  mass  matrix. 
An  eigenproblem  solution  would  then  yield  a  set  of  orthogonal  modes.  McGrew 
(Ref,  2)  proposed  using  the  Gram-Schmidt  orthogonal ization  procedure  by 
modifying  it  to  include  mass  weighting. 
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Targoff  (Ref.  3)  presented  a  procedure  which  yielded  orthogonal  modes  by 
making  minimum  changes  to  the  measured  modes.  Subsequently,  Baruch  and  Bar 
Itzhack  (Ref.  4)  demonstrated  that  Targoff's  results  could  also  be  obtained  by 
minimizing  the  mass-weighted  difference  between  the  measured  and  orthogonal ized 
modes.  The  minimization  was  performed  subject  to  mass-weighted  orthogonality 
constraints.  In  References  3  and  6,  additional  refinements  to  the  procedure 
of  Reference  4  were  introduced.  Of  particular  value  is  the  ability  to  ortho- 
gonalize  measured  modes  relative  to  each  other  and  to  another,  already  ortho¬ 
gonal  set  of  modes. 

An  attractive  feature  of  the  Reference  3  procedure,  often  referred  to  as 
the  symmetric  correction  procedure,  is  that  the  changes  made  to  the  measured 
data  are  a  minimum.  Experience  with  numerous  spacecraft  mode  sets  has  demon¬ 
strated  that  if  the  measured  modes  satisfy  certain  orthogonality  requirements 
(i.e.,  the  off-diagonal  terms  of  the  unit-normalized,  generalized  mass  matrix 
are  less  than  0.10),  then  the  differences  between  the  measured  and  orthogon- 
alized  modes  are  small.  However,  as  the  modal  contamination  increases,  the 
required  changes  obviously  become  larger. 

The  Reference  3  procedure  assumes  that  the  lack  of  orthogonality  between 
any  two  modes  should  be  corrected  by  splitting  the  error  equally  between  the 
two  modes.  This  approach  is  adequate  for  small  modal  contamination,  since  the 
changes  thus  made  are  small.  This  approach  is  also  preferable  to  other 
approaches  if  the  contamination  is  relatively  large  and  the  cause  of  the 
contamination  cannot  be  established.  Note  that  this  approach  has  been  used 
successfully  in  a  large  number  of  programs  in  which  the  measured  modes  were 
used  as  the  normal  coordinates  of  the  spacecraft  dynamic  model. 

A  large  number  of  spacecraft  mode  survey  tests  are  performed  using  multi¬ 
shaker,  sine  dwell  test  procedures  (e.g..  Ref.  7).  Frequently,  with  complex 
spacecraft,  the  number  of  closely  spaced  modes  exceeds  the  number  of  available 
shakers.  Under  these  circumstances  it  is  often  possible  to  obtain  accurate 
frequency  and  damping  measurements,  even  though  mode  shapes  of  acceptable 
quality  might  not  be  obtained. 
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If  these  mode  shapes  are  to  be  used  as  normal  coordinates,  they  will  have 
to  be  orthogonal ized.  However,  before  using  a  procedure  such  as  that 
presented  in  Reference  3,  it  would  be  advantageous  to  improve  the  mode  shape 
orthogonality  by  considering  the  degree  of  isolation  that  existed  when  each 
mode  shape  was  measured.  It  is  the  purpose  of  this  report  to  introduce  an 
identification  procedure  that  accomplishes  the  above-mentioned  objective. 
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II.  THEORETICAL  DEVELOPMENT 


It  is  reasonable  to  expect  that  more  accurate  mode  shapes  will  result  if 
the  degree  of  mode  isolation  present  during  measurement  is  used  to  improve  the 
mass-weighted  orthogonality  of  the  measured  modes  prior  to  analytical  ortho- 
gonalization.  Assume  that  for  each  measured  mode  the  natural  frequency, 
critical  damping  ratio,  and  shaker  force  levels  and  locations  are  known.  In 
addition,  the  shaker  locations  need  to  be  included  as  degrees  of  freedom  in 
the  measured  mode  vectors.  Furthermore,  it  will  be  assumed  that  any  deviation 
from  mass-weighted  orthogonality  is  due  to  modal  contamination  and  not  errors 
in  the  mass  matrix. 

The  mode  identification  procedure  will  be  derived  using  constrained  mini¬ 
mization  theory.  To  minimize  changes  to  the  measured  data,  the  difference 
between  the  measured  modes  [$m]  and  identified  modes  [$]  will  be  minimized. 
Thus,  the  error  function  e  will  be  defined  as 


e  =  ||(W]([+]  -  [0)11 


V*jk 


A®  \  ^ 

■  V ) 


(1) 


where  [W]  is  a  square,  positive  definite  weighting  matrix.  Note  that  if 
1/2 

[W]  =  [M]  ,  the  error  function  e  reduces  to  the  function  used  in  Reference  A 

and  would,  therefore,  be  consistent  with  the  orthogonalization  procedure  of 
Reference  3  (Ref.  4).  However,  for  the  derivation  presented  herein,  we  shall 
allow  [W]  to  be  any  square,  positive  definite  matrix. 


To  introduce  the  degree  of  isolation  that  existed  when  the  mode  shapes 
were  measured,  cross-orthogonality  constraints  will  be  imposed  on  the  minimi¬ 
zation  of  Eq.  (1).  These  constraints  can  be  derived  from  the  forced  response 
equation  for  a  linear,  elastic  structure.  By  taking  advantage  of  an  element- 
by-element  (scalar)  multiplication  operator  0,  described  in  Reference  8,  the 
forced  response  equation  can  be  written  as 


t4>m]  =  (<*]  ([QI  ©  ([*]T[f] )) 


(2) 
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where  the  operations  defined  in  the  equation  must  follow  the  nesting  of  the 
parentheses . 


In  Eq.  (2),  [4>]  represents  all  the  mode  shapes  of  the  system  and  thus  has 
an  infinite  number  of  columns.  The  matrix  [$]  in  Eq.  (1)  is  a  column  subset 
of  [$] .  The  matrix  [Q]  is  the  imaginary  component  of  the  system  admittance 
matrix  where 


Q.  .  = 


2C.a?. 

i  n 


ij  [1  -  a?.]  +  [2£  .a.  .  ] 2 


a.  .  =  u  ./u  .  (4) 

lj  j  ni 

The  columns  of  [F],  which  can  all  be  different,  are  the  measured  force  levels 
used  to  establish  the  corresponding  columns  of  [$m] .  Therefore,  each  column  of 
[4>]  [F]  represents  the  modal  forces  present  during  measurement  of  the  corres¬ 
ponding  columns  of  [4>m] .  Note  that  when  [$m]  is  normalized  to  yield  unity  on 
the  diagonal  of  [$m]T[M] ($m] ,  the  columns  of  [F]  need  to  be  scaled  accordingly. 

T 

Premultiplying  Eq.  (2)  by  [$]  [M]  and  taking  advantage  of  the  mass-weighted 
orthogonality  exhibited  by  normal  modes,  i.e., 


[M]  {♦>.  -  1.0  if  i  -  j 
=  0.0  if  i  4  j 


we  obtain  our  cross-orthogonality  constraints 

[4>]T(M]  [$m]  =  [fl]  e  <[4>]t[f])  (6. 

The  method  of  Lagrange  Multipliers  (Ref.  9)  will  be  used  to  incorporate, 
in  the  minimization  of  e,  the  constraints  defined  by  Eq.  (6).  We  begin  by 
establishing  the  Lagrange  function  L 
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n  m 

L  =  l  l 

i=l  k=l 


5l=l  "ij  -  V 


m 


+  Ji  i.i Xfti  l  k!i  ~  °ia  j!i  V 


(7) 


where  X^  are  Lagrange  multipliers.  Next,  we  take  the  partial  derivative  of  L 
with  respect  to  each  of  the  unknown  These  derivatives  are  set  equal 

to  zero  to  establish  a  set  of  n  x  m  equations  that  the  <J>„  must  satisfy 
for  L  to  be  a  minimum.  Expressing  these  equations  in  matrix  notation,  we 
obtain 


2[W]2([<|>]  -  [4>m]  )  +  [M]  [4>m3  [X]T  -  [F]  ( [Q]  0  [X])T  =  [0]  (8) 

Equations  (6)  and  (8)  represent  the  complete  set  of  equations  needed  to 
determine  the  unknowns  [ <J> ]  and  [X]. 

We  begin  our  solution  by  casting  Eq.  (6)  into  a  more  convenient  form.  By 
taking  advantage  of  the  element-by-element  operator  properties  (Ref.  8),  we 
can  write  Eq.  (2)  as 


I4>m]  =  l  {*M4>};[f][Q] 

j=l  J  3  J 


(9) 


where  the  elements  of  the  diagonal  matrix  [fij^  are  the  elements  of  the  jth  row 
of  [Q],  i.e.. 


Sk  =  njk  for  *  "  k 

=  0  for  SI  ^  k 


Next,  we  transpose  Eq.  (9)  and  post-multiply  by  [M]  [<ji]  to  obtain 


[*m]T[M][<|>]  =  l  [Q]  .  [F ]  T{<t»}  .{$}T[M]  [$] 

j=l  3  J  J 


(10) 
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Taking  advantage  of  Eq.  (5),  we  can  reduce  Eq.  (10)  to 

[$m]T[M]  [  <(>  ]  =  l  [Q].[F]T{4»  {I>*  (ID 

j=l  3  J  J 

where 

{I}T  =  [Ml  [♦]  (12) 

and 

iu  - 1  if  i  =  j 

=  0  if  i  *  j 

Note  that  unlike  Eq.  (9),  which  involved  all  the  normal  modes  of  the 
system,  Eq.  (11)  contains  only  those  modes  for  which  there  are  measurements. 
This  is  consistent  with  our  understanding  that  the  only  modes  that  can  be 
identified  are  those  for  which  we  have  data. 

_2 

Proceeding,  we  premultiply  Eq.  (8)  by  (1/2)[W]  and  solve  for  [<j>] 

[♦1  -  I*"1]  -  |(W]"2[M]Um][X]T  +  |[Wr2[F]([Ql  0  [\])T  (13) 

The  jth  column  of  [$]  is,  therefore, 

{♦}j  =  {  N>m]  -  |[W]~2tM][*ra][X]T  +  |[W]_2[P]([Q]  0  [X])T  }{C}.  (14) 

where 

=  1  for  i  =  j  (15) 

=  0  for  i  £  j 

Substituting  Eqs.  (13)  and  (14)  into  Eq.  (11)  yields 

(4>m]T[MHy]  -  J  [Q].[F]T[t]{C}.{i}*  -  t0]  (16) 

j=l  2  2  2 

where 

M  *  [♦“]  -  |[W]_2[M][(j)ra][X]T  +  |[W]"2[F]([n]  0  [X])T  (17) 
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Equations  (16)  and  (17)  can  now  be  used  to  solve  for  [X],  and  then  Eq.  (13) 
can  be  used  to  obtain  the  desired  modes  [$] . 


The  product  {C}^  (I}T  in  Eq.  (16)  allows  us  to  solve  for  [X]  one  column  at 
a  time.  Therefore,  we  can  write  Eq.  (16)  as  a  set  of  m  equations  of  the  form 


a$m)T  [M]  -  [fl]  j[F]T)|«|<]{C}j  =  {0}  j  =  1,  2,  ...  ,  m  (18) 

Substituting  Eq.  (17)  into  Eq.  (18)  yields 


{[G]j  +  [X]  +  [H] ^([5]  0  [X])}  {C}^  =  {0}  j  =  1,  2,  ...  ,  m  (19) 

where 

[0]  =  [Q]T 

[X]  =  [X]T 

[G] ,  =  [A] T1 [Cl  , 

[H] .  =  [A]~1[B] j 

and 

[A] .  =  '  |  {  -  [6].[F]T}[W]-2tM]t4>,n] 

[B]  j  =  \  {[*m]T[Mj[Wr2[F]  -  [n].[F]T[W]"2EF]} 

[c].  =  E«t»m3T[M]  [♦“l  -  [Q]  j[F]T[^m] 


Because  of  the  definition  of  {C}^,  the  jth  column  inside  the  braces  in 
Eq.  (18)  must  equal  zero.  All  other  columns  are  arbitrary.  Therefore,  Eq.  (18) 
reduces  to 


[g] j  +  EX]  +  Eh] j([5]  ©  [£])  =  [0] 


(20) 


Note  that  the  above  equation  can  be  used  to  obtain  a  unique  solution  for  the 
jth  column  of  EX]  only.  Therefore,  to  obtain  the  complete  matrix  [X],  Eq.  (20) 
must  be  solved  m  times,  i.e.,  j  =  1,  2,  ...  ,  m. 
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Because  of  the  element-by-element  multiplication  operator  in  Eq.  (20),  each 
column  of  [X]  can  be  obtained  independent  of  the  other  columns,  i.e., 

{X}j  =  -([I]  +  [H]j[Q]j)~1{G}j  (21) 

where 

{X}j  =  jth  column  of  [X] 

{G)j  -  jth  column  of  [G] 

ank  -  *  ajk  for l=k 

»  0  for  a/k 

Once  we  have  solved  for  the  m  columns  of  [X],  we  transpose  the  resulting 
matrix  and  substitute  it  into  Eq.  (13)  to  obtain  the  identified  modes.  If  the 
modal  contamination  is  due  only  to  the  inodes  represented  in  the  measured  set, 
orthogonal  modes  should  result.  However,  if  part  of  the  modal  contamination 
is  due  to  modes  not  in  the  measured  set,  then  the  new,  identified  modes  will 
not  be  orthogonal.  Nonetheless,  they  should  be  an  improved  representation  of 
the  true,  normal  modes  of  the  system.  In  this  case,  which  will  undoubtedly  be 
the  most  common,  any  remaining  deviation  from  orthogonality  can  be  eliminated 
using  the  procedure  of  Reference  3,  i.e., 

[$]  =  [♦][©] [&J"l/2[eiT  (22) 

where 

[$]  -  orthogonal  mode  shapes 

[0]  =  eigenvectors  of 

[(3 ]  =  diagonal  matrix  of  eigenvalues  corresponding  to  [0] 
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III.  DEMONSTRATION  OF  PROCEDURE 


The  mode  shape  identification  procedure  will  be  demonstrated  by  numerical 
simulation  of  a  test  problem.  The  procedure  will  be  used  to  identify  improved 
mode  shapes  from  contaminated  modes.  The  simulated  test  modes  will  be  obtained 
from  an  analytical  test  structure  subjected  to  multi-shaker  sine  dwell  excita¬ 
tion  at  its  natural  frequencies. 

A  schematic  representation  of  the  analytical  test  structure  is  presented 
in  Figure  1.  The  springs  represent  load  paths  and  the  squares  represent 
degrees  of  freedom.  The  diagonal  terms  of  the  diagonal  mass  matrix  and  the 
stiffness  values  of  the  load  paths  are  also  shown  in  the  figure.  The  natural 
frequencies  and  critical  damping  ratios  for  each  mode  are  presented  in  Table  1. 

"Measured"  modes  were  obtained  from  the  closed  form  solution  for  a  multi¬ 
degree  of  freedom  system  subjected  to  harmonic  excitation.  For  the  test 
cases,  force  levels  and  application  points  were  selected  to  yield  measured 
modes  that  were  contaminated  by  modes  close  in  frequency  and  by  modes  outside 
the  frequency  range  of  interest.  These  shaker  locations  and  force  levels  are 
presented  in  Table  2.  Note  that  the  excitation  levels  and  locations  for  each 
measured  mode  are  distinct.  Therefore,  the  modal  forces  present  when  each 
mode  is  measured  are  different.  The  result  of  the  orthogonality  check  of  the 
measured  modes  is  presented  in  Table  3. 

To  demonstrate  the  procedure,  new  modes  were  identified  from  the  first 
two,  the  first  three,  and  the  first  five  measured  modes.  The  identification 
was  also  performed  using  seven  measured  modes.  Table  4  compares  the  two-mode 
set  measured  mode,  mass-weighted  orthogonality  to  that  of  the  identified 
modes.  As  can  be  observed,  the  identified  modes  exhibit  considerably  better 
mass-weighted  orthogonality  than  the  measured  set.  In  Tables  5  and  6,  the 
measured  mode,  mass-weighted  orthogonalities  for  the  three-  and  five-mode 
sets,  respectively,  are  compared  to  that  of  the  identified  modes.  As  can  be 
ascertained,  in  each  case  the  identified  modes  exhibit  improved,  mass-weighted 
orthogonality.  The  seven-mode  case  yielded  nearly  exact  modes. 


15 


Figure  1.  Analytical  Test  Structure 
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Table  1.  Analytical  Test  Structure  Natural  Frequencies 
and  Nodal  Critical  Damping  Ratios 


Mode 

1 

2 

3 

4 

5 

6 

7 

8 

fn(Hz) 

4.92 

5.03 

5.26 

5.52 

5.56 

6.16 

6.25 

6.81 

0.01 

0.01 

0.015 

0.015 

0.02 

0.01 

0.02 

0.02 

Table  2. 

Shaker  Locations  and 

Force  Levels 

Shaker  Location 
(degree  of  freedom) 

Mode 

Shape 

1 

2  3  4 

5  6 

7  8 

1 

4.5  -1.0 

-  0.45 

0.45 

2 

3 

1.0  -10.0 

4 

1.0 

32.0 

7.0  28.0 

75.0 

5 

1.0 

6 

-4.5  36.0 

25.0  -4.0 

-3.0 

7 

15.0 

-15.0 

8 

1.0 
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Table  3.  Mass  Weighted  Orthogonality  of  "Measured"  Modes 


Test 

Mode 

Test 

Mode 

1  2 

3 

4 

5 

6 

7 

8 

1 

1.00  0.27 

0.02 

0.02 

0.01 

0.00 

0.01 

0.00 

2 

1.00 

-0.11 

-0.02 

0.01 

0.00 

0.00 

-0.12 

3 

1.00 

-0.03 

-0.09 

-0.05 

0.04 

-0.05 

4 

1.00 

0.55 

0.00 

0.00 

-0.07 

j 

5 

SYM. 

1.00 

0.02 

-0.01 

-0.08 

1 

6 

1.00 

0.23 

0.03 

7 

1.00 

0.10 

8 

1.00 

Table  4.  Mass  Weighted  Orthogonality  of 

First  Two  Modes 

"Measured"  Modes 

Identified  Modes 

Modes 


1 


2 


1 


2 


Table  5.  Mass  Weighted  Orthogonality  of  First  Three  Modes 


Modes 

"Measured"  Modes 

Identified 

Modes 

1 

2 

3 

1 

2 

3 

1 

1.00 

0.27 

0.02 

1.00 

0.04 

0.00 

2 

0.27 

1.00 

-0.11 

0.04 

1.00 

-0.04 

3 

0.02 

-0.11 

1.00 

0.00 

-0.04 

1.00 

Table  6. 

Mass 

Weighted  Orthogonality  of 

First 

Five  Modes 

"Measured"  Modes 

Identified 

Modes 

Modes 

1  2 

3 

4 

5 

1 

2 

3 

4 

5 

1 

1.00  0.27 

0.02 

0.02 

0.01 

1.00 

0.03 

0.00 

0.00 

0.00 

2 

1.00 

-0.11 

-0.02 

0.01 

1.00 

-0.03 

0.00 

0.01 

3 

1.00 

-0.03 

-0.09 

1.00 

0.05 

-0.03 

4 

SYM 

1.00 

0.55 

1.00 

0.31 

5 

1.00 

1.00 
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The  small,  remaining  coupling  terms  are  due  to  imperfect  identification. 
Unless  all  modes  causing  contamination  are  included  in  the  identification, 
exact  correction  cannot  be  achieved.  In  practice,  modes  outside  the  fre¬ 
quency  range  of  interest  will  not  be  measured.  However,  these  modes  will,  to 
some  degree,  contaminate  the  measured  set.  Therefore,  perfect  identification 
should  not  be  expected.  Any  remaining  deviation  from  ideal  mass-weighted 
orthogonality  can  then  be  corrected  by  a  procedure  such  as  in  Reference  3. 

In  Table  7,  the  measured  and  identified  mode  shapes  are  compared  to  the 
normal  modes  of  the  test  structure.  As  the  table  indicates,  the  identified 
modes  are  in  closer  agreement  to  the  exact,  normal  modes  than  the  original 
measured  modes.  As  expected,  the  identification  becomes  more  accurate  as  the 
number  of  modes  increases.  It  is  encouraging  to  note  that,  generally,  the 
degrees  of  freedom  in  which  the  largest  changes  occurred  were  those  that  had 
the  largest  errors.  More  importantly,  nearly  correct  values  in  the  measured 
modes  were  altered  very  little. 

For  the  example  problems  presented  herein,  the  weighting  matrix  [W]  was 
1/2 

set  equal  to  [M]  .  As  demonstrated  in  Reference  4,  this  weighting  is 

consistent  with  the  assumption  of  making  minimum  change  to  the  measured  data 
that  was  made  in  Reference  3.  The  identifications  presented  herein  were  also 
performed  with  [W]  =  [I],  i.e.,  no  weighting.  However,  although  dramatic 
improvements  were  also  obtained,  using  [Wj  =  [M]  yielded  superior  results. 
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Table  7.  Comparison  of  Mode  Shapes 


"Measured" 

Modes 

* 

2 

Identified  Modes 

3  5 

7 

Normal 

Modes 

1.03 

1.01 

1.09 

1.08 

0.96 

0.96 

0.36 

0.36 

0.35 

0.35 

0.35 

0.35 

0.32 

0.34 

0.34 

0.33 

0.32 

0.32 

{♦>1 

0.44 

0.42 

0.42 

0.42 

0.42 

0.42 

0.40 

0.38 

0.38 

0.38 

0.39 

0.39 

0.42 

0.42 

0.42 

0.41 

0.41 

0.41 

0.33 

0.31 

0.31 

0.33 

0.32 

0.32 

7.92 

8.83 

8.83 

9.06 

9.16 

9.18 

0.44 

0.23 

-0.56 

-0.41 

0.45 

0.47 

0.17 

0.10 

0.12 

0.12 

0.15 

0.16 

0.39 

0.24 

0.24 

0.24 

0.24 

0.25 

{♦}2 

0.20 

0.20 

0.21 

0.21 

0.15 

0.15 

0.36 

0.28 

0.27 

0.28 

0.25 

0.25 

0.02 

-0.07 

-0.06 

-0.06 

-0.03 

-0.03 

0.32 

0.26 

0.24 

0.23 

0.23 

0.24 

-16.56 

-18.57 

-18.65 

-18.65 

-19.02 

-19.63 

4.53 

1.67 

1.80 

1.81 

1.84 

0.57 

0.58 

0.49 

0.49 

0.50 

-0.47 

-0.46 

-0.43 

-0.43 

-0.43 

(♦>3 

0.30 

0.30 

0.31 

0.30 

0.30 

-0.27 

-0.27 

-0.27 

-0.27 

-0.28 

0.29 

0.30 

0.39 

0.39 

0.40 

-0.40 

-0.40 

-0.44 

-0.43 

-0.44 

-3.79 

-3.66 

-4.23 

-4.26 

-4.33 

-2.79 

-1.18 

-2.08 

-2.19 

-0.54 

-0.48 

-0.41 

-0.44 

-0.40 

-0.45 

-0.52 

-0.55 

{♦>4 

0.06 

-0.06 

0.02 

0.02 

-0.01 

0.04 

-0.00 

-0.00 

0.64 

0.63 

0.45 

0.47 

0.33 

0.35 

0.49 

0.52 

-3.20 

-3.17 

-2.22 

-2.34 

-3.79 

-3.79 

-2.63 

-2.74 

-0.57 

-0.57 

-0.49 

-0.51 

0.26 

0.41 

0.43 

0.45 

{♦}5 

0.05 

0.10 

0.02 

0.02 

-0.02 

-0.04 

-0.00 

-0.00 

0.67 

0.52 

0.52 

0.54 

-0.35 

-0.38 

-0.46 

-0.47 

-3.15 

-2.30 

-2.35 

-2.44 

*Number  of  modes  used  in  identification. 
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IV.  SUMMARY 


A  mode  shape  identification  procedure  to  improve  the  isolation  of  measured 
modes  has  been  introduced.  The  procedure  was  derived  by  minimizing  the  differ¬ 
ence  between  measured  and  identified  modes.  The  minimization  was  performed 
subject  to  cross-orthogonality  constraints,  which  allows  the  procedure  to  take 
into  account  the  degree  of  mode  isolation  present  during  measurement.  A  sig¬ 
nificant  feature  of  the  procedure  is  that  each  measured  mode  can  be  established 
using  different  excitation  locations  and  force  levels.  This  allows  the  proce¬ 
dure  to  improve  the  isolation  of  modes  measured  with  multi-shaker,  sine  dwell 
test  procedures.  Another  attractive  feature  is  that  changes  made  to  the  test 
data  are  minimum.  This  report  presents  the  theoretical  formulation  and 
demonstrates  the  procedure  by  an  example  problem. 
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